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ABSTRACT 

Measuring g-mode pulsations of isolated white dwarfs can reveal their interior properties to high preci- 
sion. With a spectroscopic mass of w 0.51Mq (logg = 7.82), the DAV white dwarf HS 1824+6000 is near 
the transition between carbon/oxygen core and helium core white dwarfs, motivating our photometric search 
for additional pulsations from the Palomar 60-inch telescope. We confirmed (with much greater precision) 
the three frequencies: 2.751190 ± 0.000010 mHz (363.479 sec), 3.1 16709 ± 0.000006 mHz (320.851 sec), 
3.495113 ±0.000009 mHz (286.114 sec), previously found by B. Voss and collaborators, and found an addi- 
tional pulsation at 4.443120 ± 0.000012 mHz (225.067 sec). These observed frequencies are similar to those 
found in other ZZ Ceti white dwarfs of comparable mass (e.g. \ogg < 8). We hope that future observations 
of much lower mass ZZ Ceti stars (< 0.4Mo) will reveal pulsational differences attributable to a hydrogen 
covered helium core. 

Subject headings: stars: white dwarfs — stars: oscillations — stars: individual: HS 1824+6000 



1. INTRODUCTION 

The hydrogen line ZZ Ceti variable (DAV) white dwarfs 
(WDs) occupy a discrete strip in the T^ff - \ogg plane 
known as the ZZ Ceti instability strip. Many groups 
have assess ed the location of th i s instability str i p both 
empirically (IWesemael et al.' '19911; 'Mukad am et al.l l2004bl: 
iGianninas et al.l 120051: ICa stanheira et al. 20 Q3) and theoret- 
ically (Brassard & Fontaina 1 1997c iWu & Goldreich '1999|; 
[Fontaine et al. 2003), and despite minor discrepancies it spans 
11,00 < Teft < 12,250 K for log^ « 8.0. iFontaine et"d] 
(|1982|) suggested that the instability strip is pure, mean- 
ing t hat all WDs w i thin t he strip are variable. How- 
ever, iMukadam et alJ (l2005h found numerous objects from 
the SDSS with associated low signal-to-noise spectra that 
were tenuously identified to be in the strip but did not 



TABLE 1 
Properties of HS 1824+6000 





Voss et al. 2006 


Gianninas et al. 2007 


ns (K) 

logg (dex) 
niB (mag) 


11 192 ±300" 
7.65 ±0.10" 
15.7 


11380 ±140'' 
7.82 ±0.04'' 
15.7 


Observed Frequencies (in mHz) 


Voss et al. 2006 
This Paper 


2.6 ± 0.4, 3.0 ± 0.9, 3.3 ± 0.4, 3.4 ± 0.8 
2.751 190 ± 0.000010, 3.1 16709 ± 0.000006 
3.4951 13 ± 0.000009, 4.443120 ± 0.000012 



vary to their observed d etection limits. Recent observations 
( ICastanheira et al.ll2007h have found some of these to be low 
amplitude pulsators. If the instability strip is pure, then it 
strongly implies that ZZ Ceti stars are a phase of evolution 
through which all DA WDs must evolve as they cool. 

White dwarfs less massive than « 0. 45 - QA7M(t, 
(logg fti 7.67 at Teff » 1 1, 500 dD'Cruz"etd] 



19961: lDominguezetal.1 Il999t IPieFrinferni et alJ 12004 



VandenBerg et al.1 120061: lPaneietal.1 120071) do not undergo a 



He core flash in the course of their evolution and therefore are 
left with a He core. Two modes of evolution can truncate the 
red giant branch evolution and prevent the He core flash: mass 
loss due to winds and mass loss due to binary interaction. In 
systems of high metallicity, mass loss due to stellar winds 
on the red giant branch can be signifi cant enough to lose 
the H envelope prior to the c ore flash dD'Cruzet all [19961: 
[Hansen 2005 HKilic et al.l l2007b). Binary interaction through 
a common envelope also leads to s ignificant mass loss 
dlben & Liviolll993l: [Marsh et al.lll995h . Helium is thus the 
expected core composition for WDs below « O.45-O.47M0. 
However, little direct evidence exists of the He core. Possible 
evidence woul d be the apparent over-brightness of old WDs 
dHansenl [2005[) in the star cluster NGC 6791 dBedinetalJ 
2005 ) . Though uncerta in tie s remain dDelove & BildstenI 



" Photometrically 
mined. 



determined.'' Spectroscopically deter- 



20021: iBedin et al.' '2008a"'F), recent detection of low logg 
young WDs (Kalirai et al. 2007) makes it plausible for many 
of the old WDs to be He core. A detailed asteroseismological 
study of these low-mass WDs could provide convincing 
evidence for the core composition. Observation and analysis 
of a full spectrum of the pulsation modes in a WD can 
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Fig. 1 . — The empirical ZZ Ceti instability stiip. The circles represent systems for which temporal observ ations have been perfo rmed . Filled circles indicate 
system s not observed to vary while open circles indicate systems with observed periods. These data are from'Bergeron et al.' (2004) and'Gianninas et al. |l2005l 
1200 7). The vertical crosse s represent low-mass WDs selected from the SDSS with log^ and T^ff redetermined from MMT spectra from Kilic et al. (20073). The 
diagonal crosses are from 'Kilic et al.' (2007a) except they are reanalysis of SDSS spectra and are merely candidate low-mass WDs until better spectra can be 
obtained. The square is LP 400-22 (Kawka et al. 2006; Kilic et al. 2007a). The star is HS 1824-H6000 (see Table[T] Gianninas et al. 2007). The error bar in the 
instability strip represents the typical error for measurements within the instability strip. The dashed lines are empirical fits to the instability strip as determined 
by Gianninas et al. (2007). The solid (Panel et al. 2007) and dash-dotted lAlthaus et al. 2001) lines con'espond to cooling tracks of He WDs of the labeled mass. 
The labeled masses correspond to the following model masses (in Mp,) for Panel et al. (2007,) and .Althaus et al., 1 ,200 li) respectively; 0.19 to 0.1869 and 0.196, 
0.24 to 0.2495 and 0.242, 0.40 to 0.3986 and 0.406, 0.45 to 0.4481 fPanei et al.l2007l onlv). 



produce a wealth of information about the interior structure 
of the WD. The mean period spacing of the modes, the rate of 
change in a mode's period over time, and multiplet spHtting 
of individual modes can provide information on the total 
mass, spin rate, magnetic field strengt h, mass of H envelope, 
and core composition of the W D (ICorsico & Benvenutol 
120021: ICastanheira & Kepled lIOOSi) . This has already been 
theoretically applied to dis t inguis h between C/O and O/Ne 
core WDs by lCorsico et alJ (l2004l) . The measured change in 
an observed mode period in G117-B15 has also been used 
to con strain significantly the C/O core composition of t his 
object (iKepler et alJll9M[T99ll2000tlKepler at 1012005^ 

We plot a version of the empirical instabiUty strip in Figure 
[T|. Included are not observed to vary (NOV) systems and pul- 
sating ZZ Ceti stars from th e observatio ns of Bergeron et al. 
(|2004) and Gianninas et all (120051 l2007h . Also included are 
low-mass WDs from Kilic et al. There is a notable 

absence of low-mass (log,g < 7.67) WDs within the insta- 
bility strip. There are a few possible ZZ Ceti stars of this 
low mass that are not plotted due to the absence of spectro- 
scopi cally determined logg and T^ff measurements (IVoss et al.l 
l2007h . Also shown are the He WD cooling track s of two mod- 
els for 0.19, 24, 0.40, and O.45M0 WDs (lAlthaus et all 
1200 U iPanei et aTl 120 07). The difference between these two 
models at low mass is due to the different H envelope masses. 
lAlthaus et al.l(l2001h used a IMq main sequence star and trun- 
cated its evolution up the red giant br anch at various sta ges to 
produce He WDs of varying masses. iPanei et aTl (l2007l) used 
close binary evolution expectations for main sequence stars of 



many masses to produce He WDs of varying masses. These 
different approaches cause the different remnant H envelope 
masses that yield a degeneracy in the He WD mass and its 
position in the Teff-logg plane. This degeneracy would be 
broken in the case of a ZZ Ceti He WD where the pulsation 
mode spectrum would reveal the H envelope mass. 

HS 1824H-6000 ( hereafter HS 1824 , see Table [B was ini- 
tially observed by IVoss et aH (l2006l) to exhibit pulsations. 
Their photometrically determined log^ and Tiff pla ced its 
mass at « O.4OM0 using the tables of Althaus & Be nvenutol 
(199% . This mass was well within the theoretical expected 
mass range for He core WDs making it an excellent object 
to compare and contrast its pulsation frequencies with other 
C/O core and possible He core DAVs. However, later spec- 
troscopic measurement by Gianninas et al. ( 2007) determined 
its mass to be « 0.5 IM©, beyond the expected mass range for 
He core WDs. In ^we discuss our own observations and 
differential photometry of HS 1824. In ^we apply a Lomb- 
Scargle Periodogram approach to a non-uniformly sampled 
time series in order to obtain the pulsation frequencies of HS 
1824. In ^we report the results of our observations and anal- 
ysis. 

In ^we compare all observed ZZ Ceti periods with log,g < 
8.0. It is our hope this will yield a 'zeroth-order' approach to 
He core identification in much the same way Teff and \ogg 
measurements of field WDs identify likely ZZ Ceti stars. No 
singular distinction is currently present. 
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2. OBSERVATIONS 

We observed HS 1 824 (see Table[T]for properties) on eleven 
nights from 2006 August to 2006 October using the roboti- 
cally oper ated 60-inch (1.52 m) telescope at the Palomar Ob- 
servatory dCenko et al.l2006l) . All observations had 30 second 
exposures with dead times ranging from 20-40 seconds. To 
reduce dead time, half the CCD was read out. The large vari- 
ance in the dead time was due primarily to a technical prob- 
lem within the automated observing software used to con- 
trol the telescope. The observing durations varied from 1 - 
3 hours. The primary 2048 x 2048 pixel, 11' x 11', CCD 
for the robotic Palomar 60-inch was used in all observations 
with a Gunn g filter. We chose the Gunn g filter to optimize 
the ratio of pulsation count amplitude to total stellar counts. 
It is a kno wn trend that this ratio is larg er in bluer filters such 
as Gunn g jRobinson et al.lll982lll995l) . A clear filter would 
not be optimal as it increases the total stellar counts without 
a comparable increase in pulsation count amplitude, thus re- 
ducing this important ratio. Flat fielding, bias subtraction, and 
sky subtraction were performed withi n the data pipeline of the 
Palomar 60-inch Telescope Archive dCenko et alj|200d) . The 
sky subtraction was done as an inaccurate scalar value and for 
our purposes was added back into the data and recalculated 
using standard, more accurate IRAF' tools. 

2.1. Data Reduction 



We used the IRAF package VAPHOT dDeeg & Dovld2001h 
to dynamically determine optimum aperture sizes as a func- 
tion of seeing for our photometry. Given a characteristic 
frame for each night, VAPHOT calculates an optimized aper- 
ture using a PSF for each star that maximizes the signal to 
noise within the aperture. This optimized aperture is then 
found as a function of the seeing. Then, for a time series 
of frames, VAPHOT calculates the seeing value for each indi- 
vidual frame and scales the optimized apertures accordingly. 
Finally, all aperture information is input into the standard 
IRAF task phot which calculates counts within the optimized 
aperture along with background noise counts measured in an 
annulus just beyond the optimized aperture. For each night 
of observation, 21 comparison stars were selected ranging 
g = 12-16 mag and along with the program star, counts and 
background information were extracted using the VAPHOT 
task. Additionally, the exposure start times for all frames of 
observation were converted to barycentric Julian dates. 

The uncertainty of this photometry for each aperture is 
given by; 



<^CCD = c+ni»« ( 1 + — ) (Ns+nI + Nd), 
nskv 



(1) 



( lHowell|[200 6'). where c is the number of integrated source 
counts in photons, ni,„s is the aperture area calculated by 
VAPHOT in pixels, riskx is the area of the annulus used to cal- 
culate the background information in pixels, Ns is the back- 
ground counts per pixel, Nr is the read noise of the CCD in 
counts per pixel, and No is the dark current in counts per pixel. 
We do not include the digitization error as it is significantly 
less than our value for the gain. For our observations, c was 
w 4 X 10^ counts for HS 1824 and « 10'*- 10^ counts for the 

' IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatory, which is operated by the Association 
of Universities for Researc h in Astronomy, Inc ., under contract with the Na- 
tional Science Foundation. |http://iraf.noao.edu] 



comparison stars, while A's ~ 40 - 600 counts, « 25 - 60, 
and Nd <C 1 counts for f,„, = 30 sec integrations. 

Additional uncertainty arises from atmospheric variability 
on spatial scales of the CCD field of view. Scintillation is a 
dimensionless measure of the flux variations of a source ob- 
served through a finite aperture (our telescope) due to fluc- 
tuations in the refractive index of the atmos phere caused by 
temperature changes. Young's formulation ( Young|[T967h of 
Reiger's theory of scintillation (Reiger 19631) gives Ssdnt = 
c^/3^3/2g-/,/Ao where So = 0.09 is a constant ( YoungI 

Il967h . d = 152 cm is the mirror diameter, X is the airmass, 
h = 1706 m is t he Palomar O bservatory altitude, ho = 8000 
m IS a constant (lYounglll967h . and A/ = 1 /f,„,. The formal 
photometric error for each star in each frame is. 



_2 _ 2 , 2 2 



(2) 



which determines the count level, c, at which scintillation 
noise becomes comparable to Poisson noise. This occurs at 
9 X 10^ counts at airmass 1.15 and 2 x 10^ counts at airmass 
2.0. Compared to the total formal error given by Equation 
(|2]l, scintillation accounts for 10-40% of the error depending 
upon the airmass (higher airmass account for higher percent- 
ages). Therefore, we are mostly limited by Poisson counting 
statistics, but scintillation can become significant at higher air- 
mass. 

2.2. Differential Photometry 

Since the atmosphere is constantly changing, differential 
rather than absolute photometry was used in the construc- 
tion of our light curves. We used an ensemble of compar- 
ison stars to reduce the noise level inherent in any single 
compar ison star. We used the weighting scheme detailed in 
ISokolo ski et al. (2001) inspired bv .GilUland & Browa(.1988i) . 
For our target star we define; 



x(i) = A 



Cp(i) 



, / = 1 , 



Cn{i) 



Cp(i) 



^ ,= I 'Etl^nCniiY 



(3) 
(4) 

(5) 
(6) 



where x{i) is the count ratio for the /'th image, Cp{i) and (7p{i) 
are the background-subtracted counts and uncertainty of the 
program star, c,„{i) and a,„{i) are the background-subtracted 
counts and uncertainty for the m'th comparison star in the /'th 
image, K is the number of comparison stars, and is the total 
number of frames in the light curve. The weights of the m'th 
comparison star, w,„, are the same for every image, while A is 
a normalization factor that gives x{i) meaning such that; 



Acp(0 = Cp(x(0-l), 



(7) 



where Cp is the mean background subtracted counts of the pro- 
gram star and Acp(/) is the difference in the total counts of the 
;th frame compared to the mean counts of the program star for 
that night. 
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TABLE 2 

Observation Dates and Pulsation Results for HS 1 824+6000 



Date (U 1) 


2. / (Ampl.) 


3.1 (Ampl.) 


3.5 (Ampl.) 


4.4 (Ampl.) 


Obs. Length^ 


Number ox 


Number or 


YYMMDD 


mHz (mmag) 


mHz (mmag) 


mHz (mmag) 


mHz (mmag) 


hr 


Frames" 


Comp. Stars'' 


060821 








4.45 ±0.10 (9.7) 


0.8 


56 


12 


060824 


2.68±0.16(8.1) 


3.15±0.16(9.1) 






0.9 


59 


13 


060827 


- 


3.08±0.15 (10.4) 






0.9 


59 


14 


060830 




3.14±0.13 (6.2) 




4.44±0.13 (5.9) 


1.1 


69 


11 


060903 










1.0 


68 


11 


060906 




3.13 ±0.14 (7.8) 






LO 


69 


14 


060909 




3.15±0.12(11.8) 






1.2 


66 


10 


061009 


2.74 ±0.05 (5.6) 


3.13 ±0.05 (7.5) 




4.45 ±0.05 (7.0) 


2.8 


137 


9 


061016 




3.10±0.03 (8.1) 


3.53 ±0.03 (5.3) 




3.5 


129 


9 


061019 


2.73 ±0.03 (7.4) 


3. 14 ±0.03 (7.2) 






4.0 


164 


7 


061021 




3.10±0.09 (8.9) 


3.55 ±0.09 (5.6) 




1.5 


90 


11 






Weig 


ited Average over Four Longest Nights 










2.73 ±0.03 


3.12±0.02 


3.53 ±0.03 


4.45 ±0.05 








Combined Data Set (See 




2.751190 


3.116709 


3.495113 


4.443120 










±0.000010 


±0.000006 


±0.000009 


±0.000012 









Note. — '-' : denotes the frequency was not observed to the required significance level of 90% 
" Includes only data used in analysis; excludes contaminated frames (i.e. cosmic ray in program star, clouds, poor seeing, etc.).'' See 0for 
discussion on the inclusion and exclusion criteria for comparison stars. There was a maximum 21 comparison stars possible. 



3. LOMB-SCARGLE TIMING ANALYSIS 

The robotically controlled Palomar 60-inch presents a few 
challenges for time domain observations. First is the vari- 
able dead time of 20-40 seconds after a 30 second exposure. 
Second, the automated observing program sometimes places 
a higher priority on other targets, thus placing temporal gaps 
in our time series. While data gaps can be addressed in dis- 
crete Fourier analysis, large variations in timing is a much 
more difficult problem that we address via the Lomb-Scargle 
periodogram approach. 

IScarglg (Il982l) defines a periodogram as a function of the 
angular frequency uj (in rad s"') as follows; 

/r . 



1 



J2i=iX(ti)C0S(uj(ti-T)) 



V 



X;,= iCOs2(tj(f,-T)) 



,.^,sin (a;(fi 



■r)) 



tan(2ajT) = 



X;^, sin(2tjf,) 
X;|IiCos(2Ljf,), 



(8) 



(9) 



where x{ti) is Acp(i) from our differential photometry 02.21 1 
and ti is the start tim e in seconds o f th e f th frame. Fur- 
ther considerations by IScargld (1198 2') and ' Horne & BaliunasI 
0986) showed that the probability distribution of power at 
frequency uj (with Gaussian white noise) is Prob(P(a;) > z) = 
e~" when the periodogram is normalized as 

P(io) = PALo)/a\ (10) 

where is t he total measured variance of x(t) over the entire 
time series. iHorne & BaliunasI (11986') also showed that for 
data with periodic signals, the normalization factor remains 
the total variance of the raw data with the signal present. 

The exp(-z) probability distribution quantifies the signifi- 
cance of any signal seen within the periodogram, allowing us 



to find the probability that the noise (presumed to be indepen- 
dent and normal) would, by itself, produce a power of z. This 
allows us to generate a false-alarm probability that states if 
we scan M independent frequencies then the probability that 
the intrinsic noise produces a power greater than z in any one 
of the frequency bins is 



Prob( Any Power > z) = 1 -(1 -e ■) 



■z\M 



(11) 



A periodic signal is thus significant to 90% over all M sam- 
pled frequencies if the false-alarm probability is 10%. Since 
our noise is not exactly normally distributed due to the pres- 
ence of unresolved pulsations which assure some correlation 
between data point the precise significance may be slightly 
lower than this. However, this will not have any consequence 
for our results. 

4. FINAL LIGHT CURVES AND PERIODOGRAMS: 
RESULTS 

All 21 comparison stars would not produce the most sta- 
ble comparison set. To determine the optimum ensemble of 
comparison stars for a given observing run, every compari- 
son star was compared to all oth er co mparison stars one at a 
time by calculating light curves 02.21) and periodograms (® . 
Such an analysis reveals consistent frames where a compari- 
son star has a count value much beyond the scatter of the nor- 
mal light curve. In these cases, that frame and comparison star 
were analyzed using standard IRAF tasks to determine what 
caused the contamination (e.g. cosmic ray strike, drift into 
bad pixel due to poor guiding). Almost always, these compar- 
ison stars were then excluded from the optimum comparison 
star ensemble for that night only. The periodograms found 
those comparison stars with consistent frequency content due 
to possible intrinsic variability. These comparison stars were 
also excluded. Individual frames were excluded when the pro- 
gram star was contaminated by cosmic ray strikes, or the en- 
tire frame was affected by a high background level or poor 
seeing. 

Comparison stars were also excluded when color-airmass 



PULSATIONS OFHS 1824+6000 



5 




O 



5000 10000 

Time (in sec) 



15000 



2 4 6 8 

Frequency (in mHz) 



Fig. 2. — Left Panels: Light curves (see j|2.2t plotted as count difference compared to mean counts of the program star on that night versus time. Counts refers 
to the total counts accumulated within the 30 second exposure. A least-squares fit of sinusoidal functions to the observed significant frequencies is also plotted. 
Points detached from data sets denote typical error bars. Right Panels: Lomb-Scargle Periodograms. (see Q Dashed lines denote power level required for 90% 
significance. Gray vertical Unes show locations of 2.7, 3.1, 3.5, and 4.4 mHz frequencies. 
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effects could not be adequately removed via a de-trending sec- 
ond order polynomial. This was done by comparing all com- 
parison stars to HS 1824 individually and looking for high 
levels of noise in the lowest frequency domain of the peri- 
odogram. De-trending with a second order polynomial is ac- 
ceptable in our situation as the periods of pulsation are much 
shorter than the hours time-scale it takes for changing color- 
airmass. The resulting optimum comparison star ensemble 
was then used to compute the HS 1824 differential light curve 
O2.2I 1. which was de-trended through the second order poly- 
nomial fitting, and the final light curve processed through the 
periodogram (® . 

Using Equation (fTTT i. we determined an observed power 
to be significant in any periodogram if the probability was 
greater than 90% (false-alarm probability less than 10%) over 
all sampled frequencies. A summary of all significant fre- 
quencies is in Table |2] where the frequency uncertainty re- 
ported is the separation of the frequency bins in the peri- 
odogram. Figure |2] shows the differential light curves and 
periodograms for our four longest data sets. The data from 
the four longest observations allow us to construct a weighted 
average to arrive at 2.73 ±0.03 mHz (366 sec), 3.12 ±0.02 
mHz (321 sec), 3.53 ±0.03 mHz (283 sec), and 4.45 ±0.05 
mHz (225 sec). The 2.7 and 4.4 mHz frequencies are con- 
firmed in three nights, the 3.5 mHz frequency in two nights, 
and the 3.1 mHz frequency in nine nights. Excess power is 
often observed in these frequencies on other nights although 
not to the required significance level (90%). 

Our final analysis combined all eleven nights of data into 
one data set. Sky conditions were not the same for all nights, 
so the individually reduced data as described above was used 
and then combined. Barycentric Julian dates must be used in 
this analysis as changes in the Earth's orbital position in the 
solar system can account for as much as a eight seconds per 
day change in light arrival time. This composite data set was 
spectrally analyzed using the Lomb-Scargle periodogram and 
the result is plotted in Figure [3] 

All four detected frequencies are recovered to our 90% con- 
fidence, however, the imprint of our window function makes it 
difficult to determine any gains in precision over the individ- 
ual nights. To address this concern we used a method of least- 
squares fitting of sinusoids at all four detected frequencies, al- 
lowing a single frequency to vary while fixing the remaining 
frequencies and minimizing x% ■ The old detected frequency 
was then replaced with this more accurate frequency. This 
was done for all four frequencies and repeated recursively un- 
til all four frequencies no longer changed values significantly. 
This method gives us accurate determinations of the frequen- 
cies, amplitudes, and phases of the four detected frequencies. 
To determine the precision of these new measurements we 
used a more robust X/,, minimization technique allowing all 
parameters to vary, now including the frequencies. The in- 
herent non-linearity of the fitting model requires the use of 
the Levenberg-Marquardt method which is given the accurate 
determinations of the frequencies, amplitudes, and phases of 
the four detected frequencies as a starting point. This method 
incorporates the calculation of the covariance matrix which 
in turn gives us a measure of the precision of each parame- 
ter of the best-fit model. This yielded more precise values 
of 2.751 190 ±0.000010 mHz, 3.1 16709 ± 0.000006 mHz, 
3.495 113 ±0.000009 mHz, and 4.443120 ± 0.000012 mHz, 
more than a 1000 fold increase in precision. 

This new fitted sinusoidal function was then subtracted 
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Fig. 3. — Top Panel: Lomb-Scargle Periodogram for combined data set 
(see 0. Bottom Panel: Lomb-Scargle Periodogram for de-signaled (using 
only the four significant to 90% detected frequencies) data set . For both plots 
the dashed lines denote power level required for 90% significance. 

from the data and its periodogram can be found in the lower 
panel of Figure|3] The striking features of this de-signaled pe- 
riodogram is the remainder of two signals of significant power 
near 4.44 mHz and 5.75 mHz. However, the false-alarm prob- 
ability arguments are not valid in a data set where signals have 
been removed artificially. Interestingly, the excess power near 
4.44 mHz is in a frequency bin significantly offset from our 
reported detected frequency. If we treat both of these left over 
frequencies as real pulsation frequencies and use our algo- 
rithm, we find that there may exist two more detected frequen- 
cies at 4.450643 ±0.000017 mHz and 5.755451 ±0.000018 
mHz. The existence of these frequencies is questionable be- 
cause neither frequency had enough power to reach our re- 
quired false-alarm significance level in the full data set. Ad- 
ditionally, the new 4.450643 mHz frequency is entirely lost 
within the window function around the original detected fre- 
quency. When we de-signal the entire data set with the least- 
squares fitted sinusoidal function including the six frequen- 
cies, the resulting periodogram no longer contains any fre- 
quency bins with significant power The existence of these 
two pulsation frequencies is uncertain until better data with 
higher frequency sampling and more amicable window func- 
tion can be obtained. 

5. CONCLUSIONS 

We have successfully detected four pulsation frequen- 
cies (periods), 2.751190 mHz (363.479 sec), 3.116709 
mHz (320.851 sec), 3.495113 mHz (286.114 sec), and 
4.443120 mHz (225.067 sec), in multiple observations of HS 
1824-1-6000. There are also two possible pulsation frequen- 
cies (periods) at 4.450643 mHz (224.687 sec) and 5.755451 
mHz (173.748 sec). With these periods of pulsation in HS 
1824, the question remains if it, or other low gravity systems, 
can be empirically distinguished from the normal C/0 core 
ZZ Ceti population. To answer this we compiled all known 
ZZ Ceti stars with published pulsation periods and spectro- 
scopicaUy measured gravities of logg < 8.0. This search re- 
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Fig. 4. — Spectrum of reported pulsation periods for published ZZ Ceti systems with spectroscopically determined \ogg < 8. HS 1824 is highlighted in gray 
with its four observed period locations mai'ked with four vertical lines. Some marks represent more than one (very closely spaced) observed pulsation period, 
see references for details. HL Tau 76 lists only the verified independent pulsation modes of Dolez et al. 1 200 ^. Se ve ral systems listed here a re not included in 
Figure[T]as their logg and T^ff measurements are not of sufficient precision. References: 1 - Castan heira et al.ii200d. 2 -I Mukada m et alJ2004al 3 -FVauclair et al] 
[2000, 4 - Silvotti et al. 2005. 5 - Gianninas et al. 2007, 6 - Gianninas et al. 2006, 7 - Mullallv et alj |2005l 8 - IBergeron et alj|2004 9 - IDolez et alj|200a 10 - 
iBergeron et aL1995, 11 - .Kepler et aL2005k 12 - .Mukadam et al..2002. 13 - .Mukadam et aL2003r 



suited in 30 systems including HS 1824. In Figure |4] we plot 
all reported periods for these 30 systems. Across all of these 
ZZ Ceti systems there exist many reported pulsation periods 
ranging from 100- 1400 sec. However, it is apparent that bet- 
ter than half of the reported periods reside within the range of 
150-400 sec. The four periods of HS 1824 are indistinguish- 
able from the rest of this set of ZZ Ceti stars. Further, there 
does not appear to be any distinction between the two low- 
mass (log.? < 7.67) systems (HE 0031-5525, SDSS J2135- 
0743, ICastanheira et al.ll2006b and the rest of the set. With 
this current set of data it appears that this empirical analysis 



of reported pulsation periods is not sufficient to distinguish 
a suspected He core from a normal C /O core. However, HE 
0031-5525, and SDSS J2135-0743 (ICastanheira etal.ll2()06h 
are very close to the boundary of He and C/O core WDs and 
within the errors of their logg measurements may be C/O 
cores. 

It remains uncertain as to what degree this period spec- 
trum comparative analysis can succeed. There are two pri- 
mary differences between He and C/O core WDs that affect 
g-modes: the contrast in mean molecular weights in their 
cores, and the one fewer stratified layer in a He core object. 
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G-modes penetrate deeply into the core, so that differences 
in the Brunt- Vaisala profile there (due to the mean molecu- 
lar weight; see ' Delove & Bildstenll20p2|) significantly c hange 
the resulting mode period spectrum dArras et alj | 2006l) . The 
stratified layers of material within the WD also affects how 
different pulsation modes are trapped, drive n, and excited 
dCorsico & Benvenutol20()2HArras et al.l2006h . He core WDs 
possess only two zones of He and H, while C/O core WDs 
possess the additional zone of C/O. Qualitatively, both of 
these differences would produce differences in the mode pe- 
riod spectra, and are the subject of current theoretical work we 
are pursuing. Once these full mode calculations are available, 
we can answer whether clear differences are observable. 

Most reported systems in Figure |4] were found in obser- 
vational campaigns looking only for pulsations in an effort 
to constrain the ZZ Ceti instability strip. In most cases, no 
attempt was made to distinguish observed pulsation periods 
as independent modes, as opposed to linear combinations of 
modes. This analysis was neglected in large part due to the 
lack of extensive follow up. Our observations of HS 1824 

^ http://www.physics.udel.edu/darc/wet 
^ http://www.lcogt.net 



showed most single nights of data contain the pulsations of 
one specific period and it was a rarity to find a night of data 
with multiple pulsation periods. Ideally, very long gapless 
observations on the order of several days would address these 
problems very well. These observations could be obtained 
through the use of telescope networks such as the Whole Earth 
Telescope^ as wa s done with HL Tau 76 dDolez et al.ll2006l) 
and G117-B15A dKepler et alJI 19911 [19951) and the Las Cum- 
bres Observatory Global Telescope^. We look toward future, 
more detailed observations of many low-mass and normal- 
mass ZZ Cetis to help provide a measurable distinction be- 
tween He and C/O core compositions in WDs. 
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presentation. We thank Phil Arras for useful discussion on 
pulsations in He core WDs. This work was supported by the 
National Science Foundation under grants PHY 05-5 1164 and 
AST 07-07633. 



REFERENCES 



Althaus, L. G. & Benvenuto, O. G. 1997, ApJ, 477, 313 
Althaus, L. G., Serenelli, A. M., & Benvenuto, 0. G. 2001, MNRAS, 323, 
471 

AiTas, R, Townsley, D. M., & Bildsten, L. 2006, ApJ, 643, LI 19 

Bedin, L. R., Salaris, M., Piotto, G., King, I. R., Anderson, J., Cassisi, S., & 

Momany, Y. 2005, ApJ, 624, L45 
Bedin, L. R., King, I. R., Anderson, J., Piotto, G., Salaiis, M., Cassisi, S., & 

Serenelli, A. 2008, ApJ, 678, 1279 
Bedin, L. R., Salaris, M., Piotto, G., Cassisi, S., Milone, A. P., Anderson, J., 

& King, I. R. 2008, ApJ, 679, L29 
Bergeron, P., Fontaine, G., Billeres, M., Boudreault, S., & Green, E. M. 2004, 

ApJ, 600, 404 

Bergeron, P., Wesemael, F, Lamoutague, R., Fontaine, G., Safter, R. A.. & 

Allard, N. F 1995. ApJ, 449, 258 
Brassard, P., & Fontaine, G. 1997, in White Dwarfs, ed. J. Isem, M. Hernanz, 

& E. Gracia-Berro (Dordrecht: Kluwer), 451 
Castanheira, B. G. et al. 2006, A&A, 450, 227 
Castanheira, B. G. et al. 2007, A&A, 462, 989 
Castanheira, B. G., & Kepler, S. O. 2008, MNRAS, 385, 430 
Cenko, S. B., et al. 2006, PASR 118, 1396 
Corsico, A. H. & Benvenuto, O. G. 2002, Ap&SS, 279, 281 
Corsico, A. H., Garcfa-Ben'o, E., Althaus, L. G., & Isern, J. 2004, A&A, 427, 

923 

D'Cmz, N. L., Dorman, B., Rood, R. T., & O'Connell, R. W. 1996, ApJ, 466, 
359 

Deeg, H. J. & Doyle, L. R. 2001, in: W. Bomcki, L. E. Lasher (eds.). Third 

Workshop on Photometry, NASA/CP-2000-209614, p. 85 
Deloy, C. J. & Bildsten, L.'2002, ApJ, 580, 1077 
Dolez, N. et al. 2006, A&A, 446, 237 

Dominguez, I., Chieffi, A., Limongi, M., & Straniero, O. 1999, ApJ, 524, 226 
Fontaine, G., Brassard, P, & Charpinet, S. 2003, Ap&SS, 284, 257 
Fontaine, G., McGraw, J. T., Dearborn, D. S. P., Gustafson, J., & Lacombe, 

P 1982, ApL 258,651 
Gianninas, A., Bergeron, R, & Fontaine, G. 2005, ApJ, 631, 1100 
Gianninas, A., Bergeron, R, & Fontaine, G. 2006, AJ, 132, 831 
Gianninas, A., Bergeron, R, & Fontaine, G. 2007, ASPC, 372, 577 
Gilliland, R. L. & Brown, T. M. 1988, PASR 100, 754 
Hansen, B. M. S. 2005, ApJ, 635, 522 
Home, J. H. & Baliunas, S. L. 1986, ApJ, 302, 757 

Howell, S. B. 2006, Handbook of CCD Astronomy, Cambridge University 

Press, p. 66 
Iben, I. & Livio, M. 1993, PASR 105, 1373 

Kalirai, J. S., Bergeron, P, Hansen, B. M. S., Kelson, D. D., Reitzel, D. B., 
Rich, R. M., & Richer, H. B. 2007, ApJ, 671, 748 



Kawka, K., Vennes, S., Oswalt, T. D., Smith, J. A., & Silvestri, N. M. 2006, 

ApJ, 643, L123 
Kepler, S. O. et al. 1991, ApJ, 378, L45 
Kepler, S. O. et al. 1995, Baltic Astronomy, 4, 221 

Kepler, S. O., Mukadam, A., Winget, D. E., Nather, R. E., Metcalfe, T. S., 

Reed, M. D., Kawaler, S. D., & Bradley, P A. 2000, ApJ, 534, L185 
Kepler, S. O. et al. 2005a, ApJ, 634, 1311 

Kepler, S. O., Castanheira, B. G., Saraiva, M. F. O., Nitta, A., Kleinman, S. J., 
Mullally, F, Winget, D. E., & Eisenstein, D. J. 2005, A&A, 442, 629 

Kilic, M., Prieto, C. A., Brown, W. R., & Koester, D. 2007, ApJ, 660, 1451 
Kilic, M., Stanek, K. Z., & Pinsonneault, M. H. 2007, ApJ, 671, 761 
Marsh, T. R., Dhillon, V. S., & Duck, S. R. 1995, MNRAS, 275, 828 
Mukadam, A. S. et al. 2003, ApJ, 594, 961 
Mukadam, A. S. et al. 2004a, ApJ, 607, 982 
Mukadam, A. S. et al. 2004b, ApJ, 612, 1052 

Mukadam, A. S., Kepler, S. O., Winget, D. E., & Bergeron, R 2002, ApJ, 
580, 429 

Mukadam, A. S., Winget, D. E., von Hippel. T., Montgomery. M. H., Kepler, 

S. O., & Costa, A. F M. 2005, ASPC, 334, 459 
Mullally, F et al. 2005, ApJ, 625, 966 

Panel, J. A., Althaus, L. G., Chen, X., & Han Z. 2007, MNRAS, 382, 779 
Pietrinfemi, A., Cassisi, S., Salaris, M.. & Castelli, F. 2004, ApJ, 612, 168 
Reiger, S. H. 1963, AJ, 68, 395 

Robinson, E. L., Kepler, S. O., & Nather, R. E. 1982, ApJ, 259, 219 
Robinson, E. L. et al. 1995, ApJ, 438, 908 
Scargle, J. D. 1982, ApJ, 263, 835 

Silvotti, R., Voss, B., Bruni, I., Koester, D., Reimers, D., Napiwotzki, R., & 

Homeier, D. 2005, A&A, 443, 195 
Sokoloski, J. L, Bildsten, L., & Ho, W. C. G. 2001, MNRAS, 326, 553 
VandenBerg, D. A., Bergbusch, P A., & Dowler, P D. 2006, ApJS, 162, 375 
Vauclair, G., Dolez. N., Fu, J.-N., Homeier, D., Roques, S.. Chevreton, M., & 

Koester, D. 2000, A&A, 355, 291 
Voss, B., Koester, D., 0stensen, R., Napiwotzki, R., Homeier, D., & Reimers, 

D. 2006, A&A, 450, 1061 
Voss, B., Koester, D., 0stensen, R., Napiwotzki, R., Homeier, D., & Reimers, 

D. 2007, ASPC, 372, 583 
Wesemael, F, Bergeron, P., Fontaine, g., & Lamontagne, R. 1991, in: G. 

Vauclair & E. Sion (eds.), Proc. 7th European Workshop on White Dwarfs, 

NATO ASI Ser C, 336; Dordrecht: Kluwer, p. 159 
Wu, Y, & Goldreich, P 1999, ApJ, 519, 783 
Young, A. T. 1967, AJ, 72, 747 



